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ABSTRACT 


A program was written for the IBM 360 system of the 
Naval Postgraduate School Computer Facility to apply the 
methods of singular perturbation theory to earth-moon 
trajectories. 

Verification of results against previous work was 
carried out and it was found that eaereement Coulda be 
attained only for energy parameter values of 0.707 or less. 
No solution for higher values could be found. 

The analysis of three-dimensional orbits was then con- 
ducted within this restricted range to show the merit of 


Singular perturbation theory as an initial desten toed. 
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fe INTRODUCTION 


With the advent of space travel and of complex satellite 
Systems, a revival of celestial mechanics has appeared in 
ie@eirm science and engineering. In its present form, 
commonly referred to as orbital mechanics or astrodynamics, 
i@eecombines the theories of astronomers and scientists such 
as Kepler and Newton with modern ideas and technology so 
that man may navigate his way through space. 

Orbital mechanics problems grow in complexity as the 
number of bodies under consideration grow. In fact, they 
become so complex that only the problem of two-bodies may 
be truly solved. Solutions to forms of the three-body 
problem are accomplished under certain special conditions, 
but usually numerical approximations exist for these higher 
order problems. Szebehely [Ref. 1] states, 

"Not only is the two-body problem solved -- and 

the meaning of ‘solution’ may be different for 

astronomers, engineers, and mathematicians -- but a 

general understanding exists regarding this dynamical 

system. The problem of three-bodies on the other 

hand is neither solved nor is the behavior of the 

dynamical system completely understood. 

In order to cope with the problem of space navigation, 
many methods to determine the motion of a vehicle in space 
have been attempted. The use of two-body solutions has been 
used in the method of patched conic sections. Many pertur- 


bation techniques have been applied for two, three, four, 


and more bodies. The most common method of calculating 





Pieeceureajecuories 1s probably the use of the digital 
fonouLcr. | ihe computer is used to solve the initial-or 
final-value differential equations of motion in numerical 
form and is often time consuming and expensive. 

PPeCeLomiuneomaatt1 culty in finding solutions to trajectory 
problems, Lagerstrom and Kevorkian ie References 2 and 3 
formulated a method in which they attempt to treat the 
restricted three-body problem with a technique similar to 
that used in certain fields of fluid mechanics. This is 
tnoeemethod of solution of singular perturbations by inner 
and outer expansions and represents an attempt to obtain a 
uniformly valid analytic expression for trajectory motion 
in the restricted three-body problem. The work of Lager- 
strom and Kevorkian is aimed at giving a tool to more easily 
predict ae eirect of initial conditions on the motion of 
a vehicle as it becomes influenced by a third body in its 
field of space. This would allow the determination of 
"ballpark" data for use in more detailed solution programs 
Gumoweconciderabple savings of time and expense. 

In this paper, an attempt to use the singular pertur- 
bation method of Lagerstrom and Kevorkian to study effects 
Sener la leconditions on lunar trajectories is made. A 
brief discussion of Lagerstrom and Kevorkian's theory is 
followed by some numerical applications of this theory and 


the results are discussed. 
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A. EXAMPLE INNER AND OUTER EXPANSION MATCHING PROBLEM 
Before formulating the singular perturbation solution to 

the restricted three-body problem it is useful to look at 

a simple example of matched inner and outer asymptotic 

expansions. The example is basically that of Prandtl and 

has been used by Van Dyke [4], Ames [5] and others as a 

classical singular perturbation solution. 


If an equation of the form 
Pee a, - : es - 0 


is written and a regular expansion is attempted, difficulty 
is expected since the equation obtained when P = 0 1S. Ot 
lower order than the original. 


Consider the following equation 


paiva fo = b (1) 
Hevueiniteal conditions, 


flo) = o ae yee | 


The exact solution to this equation is 
F(x:p) = (-b)[t - eacil /( “27 \ +bx (2) 


If a regular perturbation expansion is attempted one may 


write 


ane 





Pef4 phi tp tier - : 
renee ey tt 


psf es pia ie a 


substituting into 1 
Bits ets ete Tada phe eae ob 
matching powers of BR implies 


as 


so that by integration 


f$,=bx 4 (3) 


The boundary conditions cannot be satisfied unless b = l. 


If the boundary condition at x = 0 is dropped 
F.dxe)= Ci-b) + bx 
since 
Voyet=boyste ie 1 eS 
therefore 


Foye) = bx + Ci-b) (4) 


This approximates the solution quite well except in 
the region adjacent x = 0 (commonly referred to as the 


boundary layer), see Figure 1. 


a2 





Ames states that if terms are lost or boundary con- 
Givulons sGiscarded in the outer solution they must be 
included in an inner expansion which must be developed 
(i.e., this implies matching in the overlap region). 

The asymptotic matching principle states that the m-term 
inner expansion of (the n-term outer expansion) = the n-term 
outer expansion of (the m-term inner expansion). See Figure 
2 for a diagram of the calculation sequence applying this 
Peinciple. 


Returning to the present example, the first term is 
t. = Ci-bh) + be 


for the outer expansion. The first term of the inner. 
expansion is found by a "coordinate stretching." The 
dependent enn the independent variables may be stretched 
mee oer Or separavely using physical insight or trial and 
error, however, the stretching must include essential 
elements omitted in the outer expansion procedures and must 
be matchable to the outer expansion. 

In the present example one may desire to stretch only 


the independent variable such that, 


fap) = FX, 


where 


XR = sep) 


3 : 





Using this in the initial equation 





JE | SQ AF _ Si pYb me 
oe ap an bs 2 
with 
Fo)=o , FCs) = | 
where 


StpM)ydX = Ax 


The highest derivative was lost in the outer problem 
therefore a°r/ax® must be kept for the inner problem. This 
implies 
S(prfp. # 00 as are LO 
but 

SCp)] po = O as Ae OO 
implies that the solution satisfying the inner boundary con- 
dizaon iS a multiple of X and cannot be matched. Therefore, 
the only possibility is that su) //A = a constant as 


Pea 0, and this constant is assumed to be l. 


Thisce gives 
Ate d F b 
ameae + —_— =— 6 
TR ae ‘ (6) 
with 


Floy= 0, FOJpy=l , SCH) = 
In the inner variables X = X/ A and F, the following expan- 


sions are produced: 


R= ep JAFo t Arba + 


14 





JE ee AR + Abe : 
Monee x ax 


pti: oe Att A’ Fe oe 
we (AX ‘hx a ear 


which when substituted into 5 yield 


A7F J Ate Lo 
Aes 


end finally, 


and Foo) = oO 


Coed - ce * + (7) 


where the outer boundary condition has been discarded. 


To carry out the matching, 


l-term outer expansion = (1I-b)+b*¥ 
in inner variables = (fe bj+ 6px 
expanded for small. pr = (i-b)+ byX 
1 term inner expansion = (1-}) 

1-term inner expansion = AC\i- es ) 
in outer variables = A(\- Fale) 
expanded for small pw = A 

1 term outer expansion = A 


Pace litspecerm OulLer FO first term inner gives 


A= |I-5 
therefore 
mai uyGree * ) 


so that now 


= ee tay et + 


US) 





@eeiee« ~) + Terms 


Tl 
i 


Wert py. 


+) 
(i 


(\-b) + b¥ + No More TERMS 


-—> 
n 


to find the second term inner expansion 


Ate, a. 2s, 
Ae AX 


= b 


Ss. =c(e*41)+ bX 


therefore 
F2U-b)i—e 4) + plete st +h X14 --- (8) 


Since the exact solution for the outer variable is 


known, the matching from 1 to 2 may be skipped as it simply 


yrc.ds, 
[(-b = !-—b 
Going now from m= 2 and n= 2, 
2-term outer expansion = (i-b)+ by 
rewritten in inner variables = (Cl-b)+ prbX 
expanded for small pwr = Ci-b)t+t pbx 
2 term inner expansion = Ci-b) + pbX 
2-term inner expansion = C-Ya-e*) plete -)+bX | 
in outer variables = (0X) -€% + ale (CF ayabe | 
expanded for small pa = (1-6) —prC+bR4--: 
2 term outer expansion = (1-b)+bX -C 
rewrite in inner variables = C-b)tpbx —pC 
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matching these gives 
bx = bx —C 


Cc =o 


therefore 


. 
5 


This form of solution is necessary when no single 


) 


(\-b) +b % 


Gieiycisies ~ ) +}. X (9) 


esympcotic expansion is uniformly valid throughout the field 
of interest in a problem. This problem may arise through 
factors such as two characteristic lengths as those of the 
boundary layer thickness near the body compared to the chord 
Wemeun of the body. 

This long sample solution is presented above so that as 
the solution of the earth-moon problem is discussed an 
understanding of the procedures involved can be felt even 
though the actual solution is too lengthy to be presented 
in detail. The procedures are the same in the orbit problem 
as in the example just discussed but more complicated to 


Carry Out. 


B. NATURE OF THE PROBLEM 

With the method of inner and outer expansions previewed, 
atvention 1S now turned to the trajectory problem in earth- 
moon space. It is beneficial to observe first the general 
problem which must be studied, the three-body problem. 

The three-body problem is one aspect of the n-body 


problem of mechanics. In earth-moon space the problem may 


aay, 





be depicted as in Figure 3. This figure represents the 
restricted three-body problem which is the second simplest 
of the many-body problems. The restricted three-body prob- 
lem, however, is still unsolvable in a closed analytic 
manner, The only n-body problem which has been solved is 
the two-body problem. It is mrenersre of interest to 
observe the nature of the three-body problem and determine 
what information one may find from determinable quantities. 
Using the system of three masses depicted in Figure 4 
fms DPOssible to define three position vectors in space 
by Py» Pos and r3. Now, if the conventional methods of 
vector mechanics are applied to the system and the necessary 
algebraic manipulations are carried out it is possible to 
determine ten of the necessary eighteen integrals of the 
total SOTunion pouune  taree—-pody problem. 


Writing the equations of motion with reference to 





Figure 4 
mM, Y; = GmMm, if + GM, 2 is 
Yi “is 
Mls =G ee ee + GM, Ws a (10) 
¥r5 Yas 


ee os =~ 
¥3% = QMmym, [7 Gmeym,, Ws2 
ee | i—y+ aaa a | 3 ree 


31 Y32 
using the fact that 
= = < ae 
N= - Nye and an = Tye 
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and adding the three equations 


mir, Bs Ma Fe ct Mg Ts = GGnm, Te + Taira + wa, S2) 
wie Vi3 ae 


t _v a 
+ Mims tt t may, (fis) 2 m3 2 (_ 2) 
v. Yi Vas 


tt 


which yields 
Pa Gl = 0 (11) 


which is the conservation of linear momentum. 


This vector equation integrated twice then gives six of the 


ten available integrals. 


Now using the relation 


— 


YXNG HERG -Tont = HX; 
on equations 10 the following equations are obtained: 


Mm. Ch nh 2G Mee (Hy xM) 4a Mime (4 x%) 
ta ‘3 


Ma (F, xt, ) =G eae.) + cae Ct. XY; ) 
3 


| 
Ms Lisx¥,) = GM (7yxH) + OMame (¥; x72) 
7, 22 


If the vector relationship 


Ue BAO = —-; xe 


is applied and the equations are added the following equation 


is found: 


~ as MMe 
2 mi l¥o xf) = a{ To le, eA -Laxe ete ( 


wy) 





wieem Feduces to 
Smilixte) = © 

motcang that 

A (Px¥) = tx¥ 4 xF =x 
the equation 

=m. OK xe ) = Oo 
may be rewritten as 

Jmilhex% Jat = milnixh®] 
Wimeety finally yields 
S im be ae eee (12) 

wnich is conservation of angular momentum. 
This allows the determination of three more of the necessary 
integrals for the three-body problem. 

The Po owine series of algebraic equations will lead 
tomune determinacrion of the final integral for the three- 
body problem. Writing Newton's law and operating on it 
Onewoocains: 


ae 
— 


R= mt 


Ee BAe Ct. oe ) 


At 
{Fe aos VY m 
] 
Lomo = \ Fi AV; 


Using this operation on Equation 10 
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Y PoOat va e " 
Lmu, =Gumm.z ——— 44 Mm, a AY 
V2 an 


2 Yr «ay = = 
if MW2vr = G room| __ eg A ee aman) ses 
= Y,3 es Se 

zt Ti 4 


2 Fa gn Fe ~ 2 AF 
tm3vu3 = Gm; | ee AY: + mae ATs 


on ‘ Wace 


the equations are now added to find 


t 


Yc 


t a eee, ae mae 7. 
22, aA = am. ms) - At; amen ae 
tj x 


Sera \ AV, 


7 
= Ga WM. WM ; 
iy 
so that 
1 : t wm, ™, 2 
eae We ames | = ConsT (13) 


" 
This is the tenth integral which represents the total 
conservation of energy. With eighteen integrals necessary 
for solution and only ten integrals that are defineable 
further pursuit of an analytic solution of the three-body 
problem is impossible. 
If a restricted three-body system is assumed one further 


integral may be found to offer new information. This 


el 





iepral, called the Jacobi integral, does not solve the 
problem in analytic form however. 
The following two assumptions are adopted to formulate 
@ restricted three-body problem: 
1) m 


<—m- and m 


3 el 2 


2) m, and m, are in circular orbits about each other. 


With this preliminary introduction to the three-body 
and the restricted three-body problems concluded attention 
Will now be turned to a singular perturbation solution to 
the restricted three-body problem formulated by Lagerstrom 
and Kevorkian [2]. The development will be outlined in 
notation of Reference 2, however for more detailed dis- 
CUssion and mathematical details the reader is referred 


to Ref. 2. 


C. FORMULATION OF THE RESTRICTED THREE-BODY PROBLEM 

With a brief neon sdon of the general three-body 
problem completed and with an example problem of matched 
inner and outer expansions discussed, attention is now 
turned toward a discussion of Lagerstrom and Kevorkian's 
work in the application of this expansion solution to a 
restricted three-body problem in earth-moon space. 

The planar restricted three-body problem may be formu- 
lated considering the notation used in Figure 5. In this 
figure x-y is an earth centered frame of reference and x-y 


is a moon centered frame of reference. 
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Motion in the x-y system may be considered to be 
Keplerian about the earth with only small perturbation 
wemms Gue to the moon until the particle comes into the 
influence of the moon's gravity field. At this time the 
peruurbation from the moon is no longer small, hence the 
second characteristic length which results in the singular 
perturbation problem and allows the use of the boundary 
layer "stretching" type solutions. The effect of the moon 
destroys the uniform validity of the Keplerian solution 
egoureune Garth. If the mass fraction of the earth-moon 
System is used as a scale factor one may discuss the motion 
near the moon in terms of the inner or x-y stretched 
Ceoomdinate system. 

Considering the outer problem first it is possible to 


Voereven vie planar equations of motion in x-y space directly 


as 
aus cl =p*) = = pf 
(14) 
Ge tty = 8 
where 
Fo -x mo 
f = re — Fm and Y= tet tm 


If the terms on the right hand side of the previous 
equations are set to zero, one has the Keplerian or two-body 
equations near the earth. The Keplerian integrals will 


change little near the earth. These integrals are those 
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of energy, angular momentum, and position. Recalling the 
equations for the three-body problem discussed earlier one 
has, 

amc. xU.) = Const 
and 


Zz (4m omni ) = Cons T 


Rewriting these in two-dimensional coordinates in the 


notation of Lagerstrom one has 


\ ye a i 
ele at] - ISP scour = he (15) 
and 
ya i eos T= De (16) 


Two equations specifying the position of the orbit in 


the x-y plane may also be written as 


Pe = Ne X + Ci —p)4lr 
Gi) 


Ze = —Mey + Ci—-py*/r 


These are written such that if pee O the motion is 
symmetric about the x-axis. Also if q,>0 the 2 DOSEN 
of the perigee is negative for Re x 0 and always non-negative 
fOr XR = 0, 

Riecune initial Conditions are specified in terms of 
these integrals, (which will be done shortly) the motion of 
the particle in the Keplerian orbit may be determined. 

The launch is considered to be at perigee, hence 


Be. = O is given. 
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The angular momentum is defined also to be rv cos @. 
If perigee is considered $¢ = 0 and cos $¢ = 1, therefore, 
X = rv. If this is non-dimensionalized with respect to 
the angular momentum of the earth-moon system (moon's 
period) the non-dimensional angular momentum may be defined 
aS 5 | 

ve: 
= T pda lime + mn) |? (16) 
The non-dimensional energy is defined in terms of a 


non-dimensional potential energy and a non-dimensional 


kinetic energy as, 
p° = \u-py2 - U*Ph elmer mm) \ (19) 


Since Nes Xe» De determine an the particle initial con- 


ditions are specified in terms of he» Deen. and q.. 


ae 
Writing these in their non-dimensional form, 
cea A, Pe= OC, 9e7O and t= OCp”) 
one is now able to specify the initial conditions of 
velocity, altitude, and position in terms of initial 
Keplerian integrals. The following Keplerian relationships 
at perigee are known: 
—Xp=Yp= ACi-e) 
x =O 
ge = am(e -b) 


2 


1 
Replacing h, by -p ,; Panoye A and p, by 0 at x = 0 and 


e 
carrying out some algebraic manipulations the initial 
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values for the coordinates and velocities become: 


xe ERAT i= te os 


eS O (20) 
AX = O 
AX 
Ady __ z Ts “hi- i z) aaa e 
an |— pr a) 
where 
Qc: G-») /2 p* 
and 


Bs i- Leeper /a Ee | - 


It is also beneficial to transform the motion equations 
to the independent variable x rather than t. One trans- 
formation is shown below and the rest follow similar steps 


to produce the transformed equations. Writing 


- days _ aylax 
air im” kt Lax 








dt A CSCC ee 
¢. (4) = eta x ae 
Le! open 
= tl tls 
= Lay ane 
a aeee 


finally 





eee | (EY 


a 


A similar development may be used to transform a’ x/at : he» 


and L so that these equations are finally available in the 


form: 
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ey 


J! 
at = 
® 
os 
x Be 
—_—_——) 
= 
=> 
R) 
“8 
nn eee nae 
\\ 
ms 
\iaeg) 
+7 
ate 
oF 
—— 


de = t\ eae REG 


D. OUTER EXPANSION 
Attempting a solution of the equations of motion for the 


particle in earth-moon space in a power series in » one has; 
WU PY = MT 4YUOY PEP) E (22) 
thx, py = to Hpy te prtiG)yd-- . (23) 


1 
The terms RY: and ae ere wunNe \exacursoluLlOns [or the 
2 
equations of motion with the perturbation terms on the right 
hand side set to zero. In common with celestial mechanics 


mMoveation, these terms are written as, 
3/ 
iis = @* Lew —-eswew)] (2h) 


Ue 


Fai Ste (x) (25) 


where 
Ei) = swf y 52 we 
MowenGie*) = 
a= Ci=p) fee 
e= [i-zpprat/apy] ~ 


ett 





1 
If to(4») and bor) are further expanded as 


$ol4j = boo ley + prtatey t+ - - (26) 


} 


WS) = wy fey + 4p, ODS oe: (27) 
where by appropriate expansions in powers of BR one finds, 


3 


took) =| 51 w'eelx- ies uel Soe | / ey, (28) 
Lo, (x) = - bel) FL2C1-ptx)x*4 3xn* | ferexd SPE Oy (29) 
Yip = TA CK- Ce xr)? (30) 


Yt, Cx) = —AP Cex +A*)\-2 p*x)/2p Le -p*x* VI = 0) 


It is also necessary to insure that the moon and the 
particle pass within the scale distance es of each other. 
igemeiase constant expansion which includes this condition 


in the analysis is, 
T=T te Te tRT 2) 
Wnere 
aipeme= to) Gt) 
Tye = Se Ct) 


and Ty is an arbitrary term which assures that the moon may 
be allowed GOwrmeeracy With the particle as desired. 

The perturbation terms ty (x) and y,(x) may also be 
determined, the details of which are long and tedious and 
may be found on pages 689-690 of Shi and Eckstein [6]. The 


final form of these equations is: 
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x 
ox) Jey T ku — ALe)| A 


(33) 


‘) 


gine x ete @aiwlew- ealde G4) 


where 
(2) = [2 "*C)~pt2)*/2 + 2Uu-preyt- (Be) 5 p2 4) pe (35) 


and 


| ° 
Ala) = - it aaaeay | | (36) 


also f (2) and & (2) are given in terms of the moon's posi- 
trommcoordinates as, 


est Ltretey ees ae: — Cos | toole) al 


fo 2) . } t+ 2*-22% Cos [vul2) rt (aie) 
gol?) = srw Lto@ -TRI]-=@ — Sin }tou(2) —Ta| (38) 


b+ @*-22 Cos [bul2) ta] (e 


lf: t, (x) and y, (x) are now expanded asymptotically as 


x > 1 the following expansions are produced: 


eee 
t,x) = (3-27) * doy Ci-x) — Ylp) + oo (39) 
Yilx) = Cae Nes CGi-x) + $(p> + 04) (40) 
where 
1 eels ce 
Y= \\ cz) Le&da)y— ace) + me hes tn Q (41) 
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/ oe \ 
‘ig JJ etee er q.cerl per pear] i en as a ae) 


with 


Q So) Suey (43) 
and 


Ula Ty enaey (44) 


1, 
With the terms t,, t,, y,, and u“y, now available the 
4 
outer expansion for t and y is known and is ready to be 


matched near the moon. 


E. INNER EXPANSION 

As the vehicle approaches the moon it comes into a hyper- 
bolic encounter with this body. The hyperbolic encounter may 
be expressed in terms of inner variables. These inner 


variables are written as, 


7 a emcee) = Az swlt-T) J. toe Re 
ya J ye 


where pois the scale factor for the inner variables and T 





is a constant which must be determined. 


The inner variable equations of motion are written as: 


e+e ~ ese +Cor ) + wU-py Cos (45) 
5 ae Pebop a 4 Sin) FRO-PYSIHP (146) 
at a Ve 


30 





where 
Yro Kt ie 


Me = eC x EDS SST age a a 
eee), 
In the hyperbolic encounter one may express the Keplerian 
Seomuelon in parametric form as, 


a Ae Uy = 
Z = ELE —Coshryeos F E(CSUY) “Sih Sve B (47) 


Te a 
4 ean ate) = 7-6 t C(6* 41) Sr hu Cos (48) 


to = Gee eee A ) (49) 


The Keplerian integrals previously discussed expressed 


for inner variables are 


m= 2lgzy + Gar ls 


~ 2 AG dy 
cs ae ba | 

(50) 
“a 7 Ax u/s 
q = 7 2G ¥ 
rs Loe + */y 


ahece Imvegrals are related to the hyperbolic elements by, 
Gavan, & = Hiezki 
/p =-€ ne ae cos 


iOmcatirv Cut tme Maveningz of the Inner and outer expan- 
sions the values of y and € as functions of x along the 


approach asymptote must be determined. These functional 


Sl 





values may be determined by allowing u + ©. As u > © the 


following functional forms of y and t are determined to be: 





wr eins AV, Suk ~~) 

A ie oy ae ae (SE (pap) 
= = => —A — 3/2 — 2x — I 

iS = 4+ GQ hog Tat Sw Y eG) (52) 


The Us and V4 terms are the velocity components on the 


approach asymptote. 


a5 Ite 
U = zhed -canyte 


It eh o> 
Vi = ~2hQki +Ceny*~ 
j42n g* 
where 
Oi atvane = seen 
and A and B are, 
Ie PU aa and Bz -P/eh 


If the negative y-intercept is considered to be Ky» 


that is to say y = -K. for x = 0, 


al 
K, = AM - BU (53) 
U, 
and 
ie eC (54) 
Ui4 and Viv the velocity components along the departure 


hyperbola, may be written also _ 
On aA A eae 
l+2rn ¥* 


pe pert Ca eG 
o ATRL 


om 
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Bee MaArtCHING OF THE INNER AND OUTER SOLUTIONS 

The matching of the inner and outer expansions must be 
carried out at this time. Returning to the example dis- 
epesed, the principle behind this matching may be recalled. 
The only difference in the two problems is that in the 
present case, since only the first _— of the moon inner 
variable expansion is needed, it is only necessary to 
rewrite the outer expansion in inner variables and then 
Garry OuL the matching. 

For the particle near the moon the x-coordinate is em 
This occurs when t = To: If the x-coordinate is specified 
oy (OKs 


oS eS ane 


ea 


this implies that 
eo oe ae 
The phase relation was written as, 
Ta tet He Tt ye = Tor poy. + OCP (55) 


bu with t = To the equation may be written as, 


og a = pe Tin + OCp) 


so that 
ae We OC) 
with 
ee %- Fm _ v1 tA Tin 
& Vs 
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therefore 
Lu 
wap te Te 


a 


If the outer expansion romwrLttcen as 


-_ 
—_ 
_— 
— 


t= toOdt+ pLtrG) ttiGr]+ --- (56) 


y= pe Pe ON eagy = ANC) be (5) 
and replacing 
ye by = I+ pCK —Tir|2? 
f= boot ty [Ca - Te Yut +b i) + sory og pl-®) 
+Y(p)] + OCp) 
mee me Aes E= ED le ~pe 


| ie 
ba (RATE )ut tose) + Cie log Ll ttogh-a 
+Y(p)-& + 00) Cae 


also itt is recalled that 52 was written 





eat ey OS 


tne cern (1 +0273? Jog pw does not occur in the later equa- 
tion. It must therefore be cancelled by some other term. 


Since the only arbitrary term in the inner expansion for t 


is 2 , a term in t must cancel the log pp term, hence, 
_ ¥, ~ 
Ce Clot) ~ log ps + Cy (59) 
and 


=i a rp ~ 
eee toate Te ENG FOU) (60) 
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y may likewise be determined to be, 


— Me toi) + hh + BLP) - Yepy -# + O01) 


Comparison of 51 to 61 and 52 to 60 shows 


U, = U 
We =—| 
miethe equations are written such that, 


ery + Cx") 
VU, 


and 


{I 


- meee ee OCIS 
“1 U 
where 


—K= Te _to0 74+ + §(0) —¥C~) 
a. 


mee CAM, — 80, 9/ VU, 


Then it is evident that 
K, aS 
Rewriting t 
a -¥,. 2 
Z ee ctu) . loy (-%) elepeOcs ) 
where 


wb 
teeny) = hte + Ye) — Cc; 
aw) 


and 





£ = ek ti ley Tis + ocx!) 
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(61) 


(62) 


(53) 


(64) 


C65) 


(66) 


(67) 





(a 
-> ~A qe ———— 
aeons Voc 
which completes the equation matching process and determines 
the unknown constants. 


These values are used in equations 50 to obtain from 


substitution 
Qh = ItuU™ 
2 = KU 
zat 
p = Kut +Citur — (68) 


¢ = KU ~ anon) 


The term T, may also be completely determined to be 


ee Gy ek oe oy — (69) 
i = U i arv~ve 
Mies toval Or composite expansion solution is the outer 


expansion plus the inner expansion minus the common part. 


It may be written as; 


YW = poy ly) + pte) + plLQe) = nia cz) | + OCs) ( 70) 


t= tolypy + ECO +p LEG) -FRI+ og) (7) 


where r 
¥2 BE -& ECV, CosB -Uy sm 5 Y/Y, (72) 
g* - (X -GE€E CosB)/U, 


LneoecOlUglonmvOmlirsy. Order alter moon passage may be 
obtained by a similar matching procedure as that used above. 
Before leaving the analytic theory it 1s wise to note 


that this procedure of inner and outer expansions may be 
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femried Out vo higher orders of ww which increases the 


accuracy as the value of the parameter pe increases. 
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ieee Las COMPUTER PROGRAM 


The computer program used in this paper was written for 
the IBM 360 system at the Naval Postgraduate School Computer 
Paecrlity. it has been run successfully on the OS/CMS time 
share remote terminal and the standard batch processor unit. 

The program is a straight forward computation with the 
use of Scientific Subroutine Package programs from the sub- 
program library. These programs are used to evaluate the 
integrals which appear in the solution near the moon. 

MVPeCOMolcauron is Carried out in double-precision 
notation to maintain as much accuracy as possible. The 
integration routine is a double-precision sixteen point 
Gaussian quadrature from the SSP library in the Fortran IV 
language. 

The necessary inputs to the program are the non- 
dimensional boundary conditions discussed in this paper. 
From these parameters the program calculates the necessary 
constants, then the integrals of the motion equations are 
computed for x = 1 (near the moon). With these values the 
program evaluates the motion in the vicinity of the moon 
and prints the necessary parameters. Control is then 
returned to the input mode and a new series of calculations 
may be carried out. If p becomes too large the machine 


proceeds to stop and ends processing. 
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The necessary equations, a test case, and the computer 
program are contained in the Appendix. A flow chart may 
be found in Figures 17 and 18. All equations are explained 
in References 7, 8, and 9 and the same notation is used as 


much as possible. 
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iV DISCUSSION OF RESULTS 
A. VERIFICATION OF COMPUTER RESULTS AGAINST THOSE OF 

REFERENCE 7 

A computer program to solve the restricted earth-moon 
three-body problem was written. The first test of this 
program was against the data presented in Reference 7. 
Figures 6, 7, and 8 show the agreement of the present 
mMesulis With those of Ref. 7. The results for y and 6 also 
agree with Table 1 of Reference 7. As a final check 
examples one through nine of Reference 7 were tested and 
comparison of results was good. 

No data of any form for values of 9 greater than 0.707 
could be made to agree with the published data of Reference 
fio “iets inability of the program to predict accurate results 
for this case of op greater than 0.707 was unresolved. The 
probable cause of the instability of the solution is due 
to the denominator of the f(x) term rapidly vanishing near 
this point. The equation for f(x) was rewritten in 
several forms, all to no avail. The equation was also 
expanded in a Taylor's series about 0. This expansion gave 
some improvement in the prediction of y and 6, but the 
improvement was not sufficient to be used for numerical 
work. 

B. PARAMETRIC STUDY IN THREE DIMENSIONAL EARTH-MOON 

TRAJECTORIES 

The program used to predict the results of Reference / 


was modified to include three dimensional orbit consideration 
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iecevanwes of 1 = O(1). The modification still led to the 
same results for planar orbits and a study for inclined 
orbits was carried out. Cases were run to predict the 
effect of variations of energy, angular momentum, position 
of the moon and orbital inclination on the space vehicle 
as it approached the moon's position. 

The present program predicts all of the important 
parameters near the moon. The effect on changes in initial 
conditions can be observed in the changes in these near 
moon parameters. One of the most important of these 
parameters and one that requires knowledge of the other 
parameters is the parameter of closest approach to the 
moon. This parameter has been chosen to show the effects 
of varying initial conditions on the moon encounter. 

Teams @ and 10 show how the initial launch conditions 
are related to altitude and velocity at launch from perigee. 
The altitude limit on perigee shown in Figure 9 indicates 
values of lambda must be less than 3.0. Figure 10 shows 
the velocity at perigee. Since conventional launch vehicles 
cannot attain values of velocity much greater than around 
40,000 kms/hr an outside limit on velocity was taken to be 
pep OCCMms/hr. This, then implies that lambda should be 
greater than 1.0. For a given launch condition the values 
of rho and lambda can be related to the velocity and alti- 
tude at perigee. The perigee launch gives the angle of 
perigee to be W and time of launch to be t = 0. Hence, 


this defines the launch conditions completely. 
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Pigure 11 shows the effect of changes in the initial 
energy, rho, relative the earth on the distance of closest 
approach. It is apparent that as the energy is reduced 
(rho becomes greater) the distance of closest approach to 
the moon is increased. The effect of orbit inclination on 
mee distance of closest approach is Bieo eViGent. a cnls 
Figure predicts an increase in the distance as the inclina- 
tion angle of the orbit is increased. 

Figure le depicts the effects of changes in the angular 
momentum, lambda, relative to the earth on the approach 
distance to the moon. As the angular momentum is increased 
(this corresponds to an increase in velocity and/or perigee 
height) the distance of closest approach to the moon is 
dnereased. The effect of the inclination of the orbit shows 
Wp in this Figure as before to indicate greater approach 
distances for higher values Of -1nelinacion. 

Figure 13 can be used to predict two effects which occur 
due to parameter changes. The effect of qT) on the closest 
approach distance is due to the variation of the moon's 
position at launch of the vehicle. For positive values of 


T the moon and the vehicle interact very closely, whereas 


il 
for negative values of a the moon passes far from the 
vehicle. 

Tinemetrect on o Shirt On thesaistance of closest 


approach may also be seen in Figure 13. Smaller values of 


the apse angle allow for closer interaction with the moon. 


ho 





The value for planar orbits is independent of the apse angle 
and therefore remains the same. 

Figures 16, 17, and 18 show the effect of these same 
Memameters On the distance of closest approach for orbits 
@meoniy Slight inclinations. As can be seen from these 
figures the effects of changes in chee lambda, and T) are 
the same as before. AS would be expected changes in the 
apse angle and small changes in the inclination angle from 
the planar case cause no significant change in the distance 
of closest approach as predicted in the planar case. 

Figures 14 and 15 show the three-dimensional geometry 
of the earth orbit and the moon passage. The inclination 
angle, i, discussed is the angle between the earth-moon 
plane and the orbit plane. The apse angle, w, is the angle 
from the node line, 2, to the apse axis. Similar parameters 


are shown for the moon passage orbit. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The method of uniformly valid asymptotic expansions 
appears to have merit as an initial ana@lvtte tool. The 
epriity of this solution to predict values of parameters 
of earth-moon trajectories with such a small expenditure of 
computer time well justifies use of these expansion solu- 
tions. A quick parametric variational study could predict 
"pall-park" initial conditions which would greatly reduce 
wie time of total solution of a trajectory problem. 

Observing the results of the analysis of this paper it 
appears that one might conclude that the closest approach 
possible occurs when the launch energy is as high as possible, 
The angular momentum relative to the earth must be low also 
for this closest approach to occur. The moon's Position 
which is fixed by qT) and w must be such that it is just 
ahead of the particle as they come near each other. 

The only other parameter which must be considered is the 
launch inclination. The minimum value of closest approach 
is indicated to occur for planar trajectories. 

From this discussion one then concludes that minimum 
distance of closest approach occurs for a planar launch at 
the highest value of energy and lowest value of angular 
momentum possible. These are attained for near earth 


launch with high perigee velocities. 
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it is recommended that further study should be conducted 
in the vicinity of the minimum energy solution region. This 
is the region of the previous mentioned instability and is 
very important in the study of low energy transfers where 
most of todays travel takes place. 

It would also be interesting to eee. some of the 
meserictions on the analysis in order to apply the solution 
to interplanetary problems. 

Kevorkian [8] has extended the analysis to higher order 
terms and 2 comparison if this result with those higher 
Sieaer solutions would also be of interest. 

Additional work is necessary before the present analysis 
Gan be applied to a full spectrum of trajectory problems. 
The present work however indicates that useful information 
may be pened Peonevnmis form Of solucrion for the restricted 
three-body problem and that this solution has merit as a 


pevential analytical tool for initial design studies. 
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Inner Two Terms 


Inner One Term 


FIGURE |. AMES, VAN DYKE BOUNDARY LAYER PROBLEM. 
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FIGURE 2. FLOW CHART OF MATCHING PRINCIPLE. 
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FIGURE 3. EXAMPLE EARTH-MOON HYPERBOLIC ENCOUNTER 
TRAJECTORY. 
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FIGURE 4, THREE - BODY GEOMETRY 
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FIGURE 5 , EARTH-MOON COORDINATE SYSTEMS FOR THE 
RESTRICTED 3-BODY PROBLEM. 
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FIGURE 7. ANGULAR MOMENTUM RELATIVE TO THE MOON. 
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FIGURE 8. DISTANCES OF CLOSEST ENCOUNTER WITH MOON. 
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FIGURE 15. MOON PASSAGE 3-D GEOMETRY 
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LIST OF EQUATIONS USED IN THE IBM SOLUTION 


The flowchart of Figures 17 and 18 shows how the logic 
Giecone computer solution is carried out. The following 


equations are those used in the solution on the computer. 


INPUT is 
(=D _ Vp D_ 
‘ = Ye 2G(He+ hw) 
fel eee 
a4 Dal Mwt Me) Se 
SeeONoie) <= ee 
i Wet Mw 


re Phese Constant 
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EVALUATION OF TERMS IN x-y (MOON) SPACE 


HRaR eh = U4+ury/[~2 

XLBAR = RY all tA Uiturts wie +Vly-8) ~UG. +) -¥1),0)|/cosz 
Aske = O& = (lank 

PRAC= Pp ~« UE +C2Ky'? 

aeee = R= g - vcenyle 

EpsBsaes S = Jy + 2h ee 


Tavsaes & = A%*Joo(pe*er)+lU-ayulicest +a™4Ti + § 
THERAR =H = ton! (-7 Citu*y*) 


WEAR =f = Sin! (-6ttu5 */F 


OM6EG4B = AQ = Jan'(C-t/u) 
EVALUATION OF VELOCITY COMPONENTS 
Ue Ui =U 
W | = W =} 


Uil= Uy = Cet ut) “[Gtv)F*-i] cose AE) Ut") 510 Cos 
Vi =e tut) [itu a] Stu =<) (t+ v*)Cos TCOSZ 


Whe ( %s)CI+u*) SINT 


Vil 


i! 


ih 


WI 


EVALUATION OF EARTH ORBIT PROPERTIES AFTER ENCOUNTER 
Psi= Y = ten Covi tt) Un | 


HE = he,, = R4Vu - [ex 
QET = Zen = =—(Vintt)Vir + | 


Ped = Pe, = Cv +1 )U 


xcev= fe = Vat) 


65 





MEOSEST APPROACH ; 
ISAC = of = Cilzr)[ tron h) he =| 
COMPUTE PARAMETER CHECKS OF LAGERSTROM AND KEVORKIAN 
PARAL = UCLY-$) = UCCFEN -— CG!) 
PAZAZ = J + OUT, —_ Soe 


EQUATIONS USED IN THE SUBROUTINES TO OBTAIN THE MOTION 


INTEGRALS “ it ; 
W002 = t30(2) =] sin eple-°2* ~2plep*2*) ‘easy 


TooPz = t-o0(2) 3 Le/24 ape 
AC PHA2. = (2) = -[2" ~pt2)*/2 + 2a p eye 2] [alas 


= foe) = Cos[tool2) —Ts | — = LF Cos [Fu.ce) 15) 
Jit aoe 2E Cos] tool) “Te } $ = 


Sim LE0e(®) -To J _ sin[ toc) -Ta] 


Got = 4.02) = . 
fe eee aay 


] 
aa Vea Perey [aen- ater] + gra ( 42 is R 


can = S = ('S 2 tice gold cco —pta]+ tigi *@ 


66 





YSIS IS BASES 


TO SOLVE Tae 
So 


S 
AL 
AL 


a PE 


VUVOOOUOUUO UU OOOOUOU OOLOLU 


ODONOOGO0CO 
HNO FOOL Ce 
CSCOCOCCOC ee! 
ONOCOOCCOCOG 
Cocccocceo 
MNeannmoarcaeoa 
Carwrwroerartanr 
OCCOCOCOCcCO 


E 


LCULATE MOTIOM OF A BODY NEAR MOON IN TH 
BOOY PROBLEM 


ai 


TARGET 


) 


Z1$ 
ASS 
a. 
N DISTANCE 


whe IT 
AmMZSun 
—*FW CW 
qIgjorw wn 
Ona CoO 
wa Chk 
YY Orn Lt 
IJINeawrwza 
ot — Ob 
fem eI <LO. 
™—OSAe a 
aM I 
= Jet Oo | 
QF = wWw 
2xx«x oO 
— 


TRE COT20 


2Zzzrrru 
(o3] 


67 


OOCWOOCOCOCOCOCOCO 
CORAIO SG OOOH OC Ces 
AIGUEQIAQIAIRCIOIO CO ON EO SF ot St 
OQOGOC CDOCCOCCO0OO 
OCOOCCCOCOC COCO 
ocnmanrtmaenoanrrcenwo 
OArWowrarrarndtran AY 
Ge@eeGeceeceeGceaGcee “ce 


BO1V436 
800450 


a 
rs 
aD) 
‘2 ad 
- ox 
a (ar | 
<I Y 
WO Cc 
had So 
N QO 
oot ale 
oa 3 
Qa % 
<I © 
oO Q 
om © 
) AI 
oO — 
om A) { om 
= p=! -C=— 
> —~On 
= } st = oF an 
aK wd mm ON H 
+O rmx NO * 
HO 2 Gram LCoS = 
=r J—Ccae =. 
x & ant e~ x 
— pl eam i 
OW Huw O + 
Oo co | OR eae) 
~O Wrox ¢ CaN) 
be ee xYOCRH on 
~~ & ~ seO~ — 
OPW Chet oF a 
He Oo~n ge << 
Hr uUCr CW oO 
Oct crete wMn ft om 
~ TL NA G—O —O 
=O. 2. NTH 4 = 
bh J Cx OKD = o 
Or<4 WOK KS —x*«xaA 


OCC e bony tee 
MO Q2G ac Kh © 
ChkH OA YK ~ar 
CHa ~~ Ke we H VCs a 
eo a™_O~we— Cnam 

“Y QwW ZFZOocno OlC~ 
CY OOe reo e=—O CO * 
WOOrD MOM o6QN— oc t+¥% 
he eim st T Pl el ON} Da a 
Wwltromw CwZ ewe wee CN 
SUuUIO ere MOreNke Cores 
DAOC eMUrre ANsvsY OFOAKt 
COs US ZHICKCCC Cad> 


AO eOOCNnNNINnNDOW ste YO + 
AaAOrFQ0 AOorTCeoOnx © 
NO etoy —-Y OOH I!O 
wm enntetw ll ~ he ON c 
OO 2Ze2Zze ji | om p} 
SOLDNCHYXDH COniOntt-aow~ 
OA eCaDawwaneCOr ac 
SSH Ter ECO ONO HN ADO!N 
OnwacKxKCOCRKaAa Oo 
UOOere ag ei Co 


aa) 


200 





COCOTCOOCOOCG 
DH DH OHVMH FINO 
SPI DOMNINMNMOW WN 
COCOOCOCOCO 
OO CODCCOOOG 
acooercnoca ne 
(ae (o'e (oa (oG (a d(at(eg(adfad(edio. 


OCOOOCOCO 
DORHAISFIONE © 
INOWODOOOLO 
COCOOCOCO 
OCOODOCOCO 
conoec ceOrg 
(eLledled (ed (ogiad(ed(ad 


fee ee OOS Seceecoege 
td 
+ We) 
tt r.) 
¥ ~N 
# = 
* uu! 
+ eon 
t + Ne 
(o) % RrWY 
Cc + Ww > - 
¢c % <I <a 
AK + wo ae 
—— 4c Je Q. 
et + to J 
aK + 0 < 
fen" ¥ Yo 
NAO t NN CG 
Or Sg cas o 
~O * Swe = 
ed x e @ 5 Cs 
Oo~ # elix k& 
om et OO a 
= S MOay a 
Caen belied aa be — Sy 
id a) ¥- oO € o 
ag # ew> ww 
+o +: Naf x 
o~ + Oe QO &- 
oz +. YY ena 
Or DI e© OMOd 
ow" *¥ eine LUO 
Nw ¥ Ors SZ eo 
wd +- © eMOO 
—O % Ne ZONE 
my # aH «O. 
Ne + Wise m~<{ 
t~— os eM &x<O 
xO t WHeO*ee Of = 
(jac % © Seria Ot MO 
Bae 2 ONE OC OF > 
w+ + CZ ome te o 
we O Dd WU NLL oe Coe a0) 
i an) H+—~ AMCxe S>- OCF F 
CQ °o NN eTWwittao es Bea 2 
Oonw KN COe Foe NG 
e~ te o Vwa«KoC Oe o 
an HH NMMTA +O CO O 
~—O te hk f= Nek De ne 
OQ HEHE KMOM2 ee — oC ce 
Y o + eS 2K Cen 
con Ke ete AC 2O CO +O «OE 
Yw i = tWOOG. ODOCtord 
Oo - *€*F§ KOM LO CO eO ec 
we fom He mNe >We emo oe om 
~ © Fete eb KK OO cO & 


we Om EOD OOm~wrHMNOOw + wt 
t—OC" KOOCNSPOCO eO O et 
m 60 SE ANI Gee be HHO 
Wm AL ee Y © mri oO & om e(hOO>O>O 
Uw Ow th OWDKYNOUWUOO «Oo CGC st 
ASE pee he ewe Oe wwe O10 WO I 
IK AW HAWWIde elu T 
LOOGRrSHRRS NOFRE YM IAI 
OTN eH Hee SA ee ee te 
AM COHORT YMOXYR- Ara Au aAOO 
Qee SUA SB = eT BESreOOOUW 
qe <I qm 


2004 


oO 
oO 


OCOGCOON0CCE 
TOPO OMNMSUN 
rrr comodo 
OCOOOCOCOOG 
COOCc OCC cCce 
Acecrncacer 
Cearococn arg 
Sececco2doo 


(RHOAFZ)) 
0))) *DLOG( 


*€2))**( 3600/2009) )) 
D 
3 
> 
2 


oe eon? 
Owe ~-COC=— 
e+e oO 
ZOZCHa* oD 
NONnNA~DN 1 
ee eealwi=~ 
MON ONACOO 
~Q1oatrcoca 
* oS eO+we 6 
KH RO RCNMNOMA 
et ece oMkK~ + 
OF OMOM Cewe 


1.-D0,SNT5,Y5) 


ww Owe Cm ONO 0D 
ON COCAO enim 
HCH ot CONw we 
OOCHOR er+he 
CG aCw COW HONE 
OxnAO+OC+~O0O0T7 

- OO FT ot NM 
ai I> YAKIMA 
AWM LT OY eee ew 
ANAnIMu Ov 


OSOCUCOK OD 
ai <a 

© 

© 

© 

t 


68 


ORBOO0860 


01 


© 


WRITE(6+104) CF11,CG11 


COO 
Ort 
One 
Coe 
OOO 
fo myeazea) 
(od [ed (oe 
Occ 
+. 
om + 
NO 
* 
%O 
mo< 
Prt, 
<q w 
+ 
«x > 
— | 
+m 
Cr 
O— 
a 
Coe 
~~ | 
be ae | 
~~ 
om UL 
NO 
Cd 
Sad 
=) =) 
> ++ 
a) (o) 
Zz CoO 
ec oO 
ot @ 
ox we \S 
Oo = 
front Cm 
- Cas 
© o~ NF 
2z2OQ qt OF 
OW oN WOm 
wNK SOx 
On Onda 
CLe=D NO]! 
Cw + we COU 
#O ZF enn 
-EO me IO 
OAD o 


a¢+xCO~~™ 
mMO~ OH SCEZNKOC 
LON eI #OOUM 
cOCxXker Het ol Y 
WHO Ute Dee Kx 
Re 9pCOOCOKNw~e+ 
Di Aes +O Wu D> 
Ad We SCeOroMnaeydo il 
Saevr aati <c~a& 
OmaAaAWUMr c— OF <T 
OF MWON 2 SMa JTDMH 
<I O Kw kK +O 
etN ot 


OOF OOUO 


ORBOO0940 


aml andl © & 
Na@W 
e-mIMN 
KOO 
we. 
=+# O 
<= Gh 
lene 
< m eH 


OMI hk awit 


+ Ao 





teZa cOnNne 
CeTA<VNCWORe 
ICre rN oLMy 
MnaTITOre ae 


Om oo mn) 
INO an Sr 
Oo (onl @) © 
Oem enim = 
Oo Oo (ep) 
cm [oaye sn) (oa) 
(oaloe [omied (ad 
OO Oo © 

@ 

(S 

LL 

= 

oO 

e 

ce 

<a 

@m 

- «I 

© ke 

[ aa) 

r.) or 

Cal) <I 

~— (sa) 

O°r Wu 

Cat ag 

em -- 

jm - 

me OK a 

~—_ <I 

KW m 

co WY 

<IU és 

~ ma tL! 
em rn a iI > 

ENDND — OO 
CFZEYU O G O 
Ie xx © -_ 2 
aD~ ¥ [sa fe5 Cat 

At Sm NSN 
—xKOQOCD Gf GQ = 
ao~CywmNe WM A 2D 

Ist «Om Coal 
Orato FT LY WN 
Tdeedeis 2 UI 
$+ ORK GSi eo Om 
OLeCOodaem a & 
OF CGOton ~ 
00M eDOt+tO X& YO 
NOONCHO <<" OC 
“ort chk c« fC 
eA Ort+O WY WwW 
C+caO~—-—OD «KX SO 

CCM amOer o 
VWOIUM+O~ aw © 
Qexew— 6K <= 2 
NA hE ONe GD = 
De~ FOND TFT = 
(2) 
O 
za 
= 


~ 


<x tmae Wow 
NCS ke Owe bE 
CdedadrdawdgD 


ADMONaAOrrOADIOZA cOlW wt 


R 


(a) OOCOOGTOO0O 
oO OM OO Oren + 
oO meter OUOIAIOIOI 
wi ated edited d aie 
(o) GoOCGeCeGeace 
a macnn nnwMecs 
(eZ (ed (ag (ee (oc (od (oA lad (adlo" 
© QOCeaeece es 
<I 
(Va) 
Y 
<I 
OF 
za 
OQ 
(2 
= 
ew 
Ui 
= 
Ww 
<I 
m i 
_— om — a 
ar & - 
qa < (ad 
cm uw <I 
dt WW 
<x OO 
sr wor a oO 
OM be 
* OO w 
~— OOo > Ww 
a QC um > 
we Ht eH — 
OF ~— << fe 
eS (OCay <I 
wit ~—OOC WwW = 
wo Aawui wf lu 
>O ASS [ad 
6. 2 2 
Lin De~ OC Pad 
O~ ¢€2N O O 
-—Or-O = — 
MCOCNL ke 
—-CO ea am oOo 
ZM ets *# W OF 
WOr~OG Kf oO 
ZH ewe URNS e 
Omnmerr Ie OO ° 
Aa-—Ne>> DWMILO© mit 
ZSNENIT t+ Sxnat CO 


Ode awu'OoO emer | 


at Ob te OO 
bg OL eC? 
PH TVW HY 


Ort o 
m= COLE 
mee 


De jue or em 
~ bh & = SO— 
—~+t¢hee te DH UWS oe 
OK CHF COM ORs eCOnO eX Oo ome oO 

CAME OM-C oKDO eD> HSHEMOOMWOL 


N 8) 

© 

N 
et ~ -~ 
~ ui OO -— 
e eo oO _ 
* it en 
t — wt 
+ m | 
+ LL —_ = 
* ee a oe) 
& GN iW) 
” ~ + a 
tL Q + 
- -<t *¥ 
ean ~ ioe — 
Oe Oo a 
CoM e -M Ww 
ah QL ma FT 
tus as £ + 
>em— UW OY CSC 
on «© Cag) 
Oed HH @ e 
ZAIN] CG) Vee BE 
met mm OH + 
Qwest Ww © O 
Cetthr © 4 
ke le meme ic 
Demi ennu ex 
OZ eCO e+ ~ 
emetgO0O kk 
Qe tr 0eO x 
mZOrzmeN o @ 
mat elUOriet WY 
=x FOAaA cU~ OO 
Ch ens eo ] 
>) 
Q 
© 
ote! 
—_ 
* 
_ 
f=) 
tu 
ae 
sk 


enh > >*e We HO oF 


OCOoOCO0O00 
COOMA SM + 
MIBLO CCE EO OD 
edt ett ete 
Cccococ 
eSfaGo@ (oai(oatab) (on, 
Orerarradr 
OCOOcoOCCe 


a nwa 
<WC~- 
ellie + 


ZOICS Sete NOH e+ KO ere & ehe WO 
KHNKMNC CHK + CU MEOOCm & &MO pea 
CAO ~~ ew Oh eH tH OODKON—OW oF 6 


OWFtunrzteZeCm erZ>o> > eo 


frm ew om ee CD es 


—C 


Ooco 
Or ao: 
cm rman 


ORBO1420 


PARAMETERS™= >?) 


-—O+x* 
<I + 
SON] ool Ow 
caga 
= YuUiel 
aI<aO 
ao Ik 
—x«K aed Y 
mw ests Ww 
Set I ASTM 
qOonwvraO 
COUOr fou 
CleHa~wO 


*2)/2.D0 


oie! KUDNOtT HD OD 
> Ci + peer et ORIN FOr He Seed Niet t+ De 


DNNNOWY RD CK FWO 4/100 Men KN HOLA O © ee 
COX OCOCMIwenZ aI owemHOOOW eMOCnmO# CKHOUTOOO oe 
OF ORY HHS AOD me MTe VN eTWwOVvYownM~AMIIT™ 

CO W eD eAQQke THT DSS emt tora excl CO eh SY ernie 


CO tt Ne ee tl 
CDR O~w mem Se TOW He 
a cOXHONOD>S>—ONs om 


mh OS Re IH HH ODK OVE 


wr De WR ee OM SH I mer 


weWtoe PRTC HH NY CWS 
ie ad oO al od ol oo on Mla Mel oo) oO <0 ae 


OAS Lh et BoB tet et te a Le a LL ot dt te tt OLN OL bem Dt se OF Cr OL OC OL 
OO FW TO RCO et te te rr TY UID IL a oO 
Wert & +h SE eOD>S Sia L ACATuU ee 


ei 


69 


eit 


3900 


u sOoCCl 
<4 


3001 


SZCaatrwweritdtearnotwr 


Zaa ZuUUU 


WC 
Ocwn 
COCe 
om 
(aa) 





Li} 
oe ta 
" We) 
reo 
<aIN CY 
OOH 
<I «all 
Nhe 
@ ot jl 
ew m 
Ce <i 
un 4] © 
be ct Us 
ote 
Oe O 
eoVWe 
Oa & 
el © 
uw Ww 
ee 
Y eS 
rewoOo 
<IAtMN o 
creo 
a or 
06 UL 
oe oe 
On Il 
Ore 
COLL <t 
hte @© 
(| > 
OM <I 
oat ke 
(aN ie a 
aCe 
Ww o 
ee © 
Wo «#9 
aor 
<aIm e~ 


ODOC Gee 
NE DOO OM 
tty NN 
te ltetediet 
OOCOXCOO 
ma otda: od 
Carvrarr 
eee ee 


/ 

Y 

°o 
T=*E1206,T50% 


) 


Whe & Tun 
© J I Oriek 
We <tett-<T us e 
aHrO elen 
wWwoOoowos il « 
QUWkr— eLiIDoO 
he Weds c 
oe 7 ecto NA 
CLO) ee et 
Oo Mel eu 
Wik IlooOe 
-- I HN UW It 
eh Ot of ei 
(0 @) oO e-- 
en Qa oO 
NNN of 
MmNLLUOHA & 
tse © WHO 
eN ff ee WIN 
EAU ile & 
=nardWil e 


MeO ecI~NCke OC 


xr * of 

CA o 
oe erin 
oO) LW ew 


e?LU fle 
ke o Il 
Com] a0 4 


OnO ee o 


NF ONAN 


THE INITIAL PARAMETERS ARE 


e eA! el & « 
Orme UW © &Ort eOer OO 


COrtUWwoni il 


et, 


CMC CALL FOR THIS CALCULATION 9>55>3>>555>5>>>°) 


Nee oF ANE ome 


cal en Le oC 


OY Mat © “CORrOCwcng 
Teqwoxx Trews + TOT 
eA OO Tete et Oe etc] fl rd oi 


ww Oh be weer wm yf me et Se OO Om 


- er EUDE SRS E 


dee Odddide e coe Tdtiagoiid 
SzSeeS eVSeSTREZTA 
CCORFAYYOOCOZHORY BMAX OC 
Cetet eOCORRHON COVTCACr2Z 
UrRNULI eR ire USE OBuU NW 


= © ere 


aco 


fol 


aa <M 
| 
(XS) 
HOM — 
oo © 
o~trn aa 


100 


Q95 


COOOCOCOCOCOCOOCOCOCCAOCOOCOCOSCCCOCCO 
DANII FIVEE OOCHAMOST NOE MOAOMNUMFMOOR MOO 
NODOODODOOOCODR ER EEERE RPE OCA MEeMMO AOS 
wh dt mt med tat mi ed pd tt mt pd at gd et dt mt ed td mt pt et 
SeeGeGoeeeceeceoeqoGgeGeeceeeeoeeec ec GGG 
LCOCGeoerawMmerOMc anne nwiacenenmad moa wa 
(0 ia Mle Mie Mia Ain Ale Mle lo Me Ae Me Mt Ale Ae Me Aa Ms Me Mle als 8 Me eae Ae! Mame ae eae Ss 
Se@CeeeceeceorveocooucecooccocoCcCoccco 


re 
e 
r) a 
e © 
(-) = pe -~ Lis 
) > oa © < cw 
° «x oc W oe ==) 
° et ee ey k= - 
© © Careless W <I 
° con >>Ww << LL cw 
e rw = Cc 
te) al (Lith mae <I 
© bd trOC 4 <I = 
r) Zan << od @ 
© > wt) a) 
© a) m © 8) VY) 
6 oc oan WiwOoO WwW fag Y 
© Ww PB Irmw ef a i 
r) UG: ee Y 2 Oca << 
C) ~< = Q (VU) =>} ies) 
c a5) ULL mY) 
— e ox a Cee rzZ Lom a 
- r,) uJ <I at al) nD) es za 
©) ry > - OCm we CO _ 
ra r) OQ Ww) Ze ay (A) UJ) © 
= ° DTG = oa ow 
a. © Q = Oleaeasy (S) .) } 
a r-) tu << mc Ww GAZ Oo 
© ° = z al} “(ed ==> = 
b- ° = (og Ard a Itt 
—”Y e = wu) lula we © (ve 
© —"Y Fe ZMmAm wu C2Y & 
Ce) e ~< Gant Ss) Oou: 
Q e iw are amt oe) ts: 
r) ° ~ <=> AWD = 
os ° x za a2zlWw oO Gre~ <= 
° ~~ ~<_ COO Ca DUOwWw lw 
za e and > rate Ww = 
<I e O CY) NAMWNATWWO (0 Ml eed 
ae ° Ww UU! One = zo. ~ 
—_ e ~ Or Wool = [es a 
.e all Che bk LLU OS me UO 
Ce e «I DD WAYY dw bk uJ 
WwW e ow Da? pzeL(eu(Uled =a} wmOQ Zz 
- C3) oO <w <1 ZW PAN UU (Ss 
<I C) tu oY wWwLhI CW Ee) (Gels Q 
Wi e oo ool [<n ee Pl lo won 
[aa e © z Kr %AarmmOm ao oad VY) 
(&) o et b— ~~) SOW ou (Qype — 
e © ici @) Gas 25 Zae 
Oo e CG Se) Ne) | BY and a) oo —aIzD> PE 
2g ¢ © t rte eaee (Ss 
or e 3) COU | } MW — 
e tw a Cre O Wp Ud om 
~e oe = Cus ZzxKxa <I 
ma Cc - (a) = & mm LY =) 
Oo e k&e UK) wt a Wt he fF at, 
Oz oD WN ae me - Ve Sn Cet 
ort oO @ (fle) Wyxsier C4 fa) cag) (ofa) ee 
© ft AkF COO OU +». We Ake SE Tu 
Ce ey =) Ue Ww Li =] tu 
oY a =) (om) ae Y) = 
ry 
e 
r.) 


oO 
oO 
CAO 
Oo 
omey OV OU OO UO UO OOOO UOOOOUOQOUOOUKOO 


(8 





OOoQaqgee 
MMO FO OF © 
OOoOrTonnro 
ee tended on 
OOTOOMOO 
canmodoce 
Crarrtraenw 
OooccoGm 


at 
AND 337-3400 


RALS» 


G 
i 


& INTE 
LOO= it 


ON O 
PP, 


I 
’ 


EO 


A 
6 


Hor 


OLYNOMTALS UP TO DEGREE 


Aww & 
wiz 
MMOD 


TE 
D 


ellie 
—ul--O 
aod 
Cc2Z2=™ 
Witte 
eo <a’ 
ZwCco 
mm Lt Ox > 
wd. 
reas 
Oo <aIw 
mer = 
ITO e 
su> s+ 


~ ec 
a>>W) 
cool ed CYS od 
Dk ew 
=O e= 
(oa GL) 
Ox ect 
Uwwos 


e2e00e20e20€¢2020002090000000909000920000000000000090000000099020000090820089000900 


OOOO OO YO 


SUBROUTINE DOQG1L6(XL»XU,FCT+Y) 


GieQG eC GeeeaG eG OoageceaqgdaGeeaceac@ 
AHOMNMSMNDOE TAO tNQ\OaASOOP CA ORMNM 
BDRADDANDDCOC CM AMA HM see HANOAICS 
HAAICUAIAIAICIOAICUAICICIAIAQIOAIAQ HIOIOGIOU II CUR GIA 
GSOGeGeceGoeaoGoeoqoeeoeececeocecGeoe 
ro Gl ° OF OS, ao Mle Die Ge Steal eS] o Gee Ole Gi aGl? Ml oAyeale Ryser Or sAl oats AloSree) 
CaMrrnnwrwarrrrrtarrtadcntortantanrad 
SCG eeceoee Cece eoeeeGeqccece 


*B 
Date CECTCA+C)4EC1(CA-C)) 


SDaTHUECTONECTSECTIA=C) } 
BN-1*(FECT(AtC)+FCTCA-C)) 
D-IUECTUARO} ECT (N=0)) 
Dt ECT(A+C) #FCTIA=C 14 
ZDSD-1e (FCT(ASC) #ECT(A-C))) 


ania mM 

¥OY - a 

ONO OOaC 
CCAaArOrCRrROaaTM 
MmPrOOMKr MAOKr OaNOr OW 
AFOIMOCONODTOAY ANE 
TOONMNFaHODOANOMOMAOKA 
NE eH OANOF CUNT AS CON 
Or ODOATAINME MOO ON RW 
MDDOOMAT™mVAHOOOOAIDO 
TINMADOMmOOCYNOTNOM 
TOD AANA COO mh ein 
MM AON VAITDOMMNOAN 
OCMEONOTAACAKE th 
TFANMWOOD tar Mor ANS 
ONE NAM NAB OOOR HOWNTM 
OQDODWAHAHNOMMEONOMO 6 


co 
Ye 
© 


DoT CF CT CASEC J SECECA-C)) 
46D-1* (FCT (AtC)+FCT(A-C) ) 


DO*B 


NO*¥B 


°(XU+XL) 


Qi Aer aM ef eC eNW est oF ~ XX 
WMOtaeteT ¢M4¢ Ot Ni ¢e¢ FH D 
ex € © O> @> o> O> O> o> 0TkO 
ot We tt th ae 
IMO>O> O> OF OF OF U>O> aw 


DOUBLE PRECISION XL eXUsY,AeB,C,FCT 


OU oO 


jer 


CoOoQoeoeceaeecqge 
TIDVOMH DO CKetlimesw 
QINIAINOLEAU EE OHM MCE 
CIAITAICICIQIQ\ QUALI AIA 
Seecocoaoeocoe 
Noaeeoroijcanrceae 
Ceraertrwrorrarera 
C@oececeeceGceGceG 


~ 
~_ 
Nia 
~ 
he Y 
[oa am} 
Cr 
om NO 
NOt 
- ew 
Io+~ 
(erp Pag 
- Con~ 
ANON 
aoe (2) CX 
NO ea 
¥# CN 
ke ont 
NO = 
+ HNO 
= se 
NO-m 
HONS 
+ eF © 
CONF OQ 
I~O o 
cH ON 
— = CY = 
(mws 
Ne 1 O 
“—# CO 
t-OO o 
er «eM cOo 
COO wt mb CO 
Dre ser er OE er eer 
Qw e- 1 O~KN AH 
ete CNOw~ 
On-Om c=—O+ 
THNH]HVUnNOe 
Noe WeROnmw~w] Ir 
1OoO #NF NRK Ka 
OF OF K¥ XCKL 
eo OF HK Cro. 
pErst OCNINICGMNOAO 
(OO Aww H OTOT 
Tou wHEH HE OC! 
awOM Zeer 1I NA 
NOT sNGOTINN 
—%+@&% NRNA¥ OOT 
aAatN MK O~WKF OTCO 
eau acs ICr ci 
zw OF-OTt~rd 
NLO ~x~orne~ 
© ~~ e- OF ¥ 


*% RHO* 


DO 
92- 


22-2500 7 -“DCOs (1e 


PTO) 
10D0-Z))) 


+(7 
=CA 
Q* ( 


7 $) 
ALPHA1+Q,TOOP1,TOO1] »YL201 »CAPT12 CAPT USN 
BT@) 
OZ 
/{ 


aka | Lao e 
Ce™ NU NO OD 
Ls VD) mee BS me et 


Ree OH FAY KO || 
Odes +Na Tee Illeana 
Z2ZoszZ2ocanenik 
DTC WOMmUAONCIZ 
the Zk Oc www OW 

aqIqm a 





OO 
or 
tO 
NN 
CoO 
mca 
(ee (oe 
OO 


RETURN 
END 


DOODDOVOMOOCOOO 
MAOMAAIMSINOM™ OOO 
SI FMIMWINU MINUTO 0 
MURIBGAIQUNI CU RICE AI AUC AIT 
MOOQOCOOCC.OCOOOCCO 
Mnrcorwmnmaca cCoOomowo 
Krannrrrrrrrrace 
@eecececoc ee co 


‘ DO * RHO* 
-CAPTO))** 


) 6 
oMO))) 
TOOL—C 


0 
0 
# (10D0/(Q%(1,00-2))) 


)) 
/2 
S ( 


C100 
# OW 
¥ 0S 


eer tee PS! 


be 


(les DO-— (RHO MS? 14773) 
if 
0 
+ 
) 
i 


KNOCOW 
~"OMmOK FEO 
HCNCHA~Y— 
Q 00 eaNatz+ 
ON OAK Onn 
O~O~ F 1 VO 
a~wO RFNRFZANCO 
~N™N =O 
O~ DN HCN 
Crary ema 
cH opt Own O 


=i, =r e) 
CAPTO,ALPHAL,Q,TOOP1L, TOOL,» Y1201_ CAPT12,CAPT ,UyN 


A-H, 9-2 
eCAPTO, 


T3(2 
AL*8 
F/RH 


awe HCHO TeZK-O 
OCO ~N~~ CeO 
(an) {Ot —WO ¥ 
zea ic oN 
CRN fH i Ne 
—~UZ I~ S# NTS 


RP eCReRtKRN FEF HO 

OWS +¢+<aacdnyry Cm>D 

ZaFZeOrFOCCN: oF§ KG 

Sa eocynwewerean2uizZ 

eS 
<a 


‘fe 





REFERENCES 


Sze vce t wm Ge. Uneory vol Orbits, p. v, Academic 
Prego. 2967. 


Lagerstrom, P. and Kovorkian, J., "Earth-to-Moon 
Trajectories in the Restricted Three-Body Problem," 


Journal :de Mecanique, V. 2, p. 189-218, June 1963. 


Lagerstrom, P. and Kevorkian, J., "“Earth-to-Moon 
Trajectories with Minimal Energy," Journal de 


Macanique, V. 2, Dec. 1963. 


Vane Dykes se berpurbation Methods in Fluid Mechanics, 
Academic Press, 1964. 


Mies. aN NOM inear Ordinary Ditferential Equations 


imleeieperpshrOceeses, p. 1f/2=-1/9, Academic Press, 
IES Ver oie 


Sie eehodesecstein. Mo. (Uniformly Valid Asymptotic 
Solutions of Non-Planar Earth-to-Moon Trajectories 
in the Restricted Four-Body Problem," The Astronomical 
Jourdain vee. pe O05-/01, Aug. 1967. 


Lagerstrom, P., and Kevorkian, J., Numerical Aspects 
or Unittormiy Valid Asymptouic Approximations for a 
Class of Trajectories in the Restricted Three-Body 
Problem," Progress in Astronautics and Aeronautics, 
V. 14, Academic Press, 1964. 


Kevorkian, J. and Brachet, G., "Numerical Analysis of 
the Asymptotic Solution for Earth-to-Moon Trajec- 
tories," AIAA J., V. 7, p. 885-889, May 1969. 


Kevorkian, J. and Lagerstrom, P., "Non-Planar Earth- 


to-Moon Trajectories in the Restricted Three-Body 
Proulene sili Jl. Ved op, LAO=152. Jan 1966. 


is 





PTT TAL SDESTRIBUTION List 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 


EGcOn. 0 mamcolinns, Code 57Co 
Department of Aeronautics 
Naval Postgraduate School 
Monterey, California 93940 


Chairman, Department of Aeronautics 
Naval Postgraduate School 
Monterey, California 93940 


LT(jg) Vernon Curtis Gordon, USN 


Route 7, Box 230 
Athens, Alabama 35611 


74 


No. Copies 





UNCLASSIFIED 


Security Classife. ation 


DOCUMENT CONTROL DATA-R&D 


(Security classification of tifle, body of abstract and indexing annotation must be enfered when the overall report Is classified) 
2@. REPORT SECURITY CLASSIFICATION 
Unclassified 

2b. GROUP 


SINGULAR PERTURBATION PROBLEMS IN EARTH-MOON SPACE 


1. ORIGINATING ACTIVITY (Corporate author) 














Naval Postgraduate School 
Monterey, California 93940 


3 REPORT TITLE 






4. DESCRIPTIVE NOTES (Type of report and, inclusive dates) 


Master's Thesis; June 1970 


S$. AUTHOR(S) (First name, middle initial, last name) 


Vernon Curtis Gordon, Lieutenant (junior grade), United States Navy 


6. REPORT OATE 7a. TOTAL NO. OF PAGES 7b. NO. OF REFS 
June 1970 > 9 


8a. CONTRACT OR GRANT NO. 9a. ORIGINATOR’S REPORT NUMBER(S) 
b. PROJECT NO. 


Cc: 9b. OTHER REPORT NOIS) (Any other numbers that may be assigned 
this report!) 
d. 


10. DISTRIBUTION STATEMENT 
This document has been approved for public release and sale; 
Peomerstribution is Wnlimited. 


11. SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 


Naval Postgraduate School 


Monterey, California 93940 





13. ABSTRACT 





A program was written for the IBM 360 system of the 
Naval Postgraduate School Computer Facility to apply the 
methods of singular perturbation theory to earth-moon 
trajectories. 

Verification of results against previous work was 
carried out and it was found that agreement could be at- 
tained only for energy parameter values of 0.707 or less. 
No solution for higher values could be found. 

The analysis of three dimensional orbits was then 
conducted within this restricted range to show the merit 
of singular perturbation theory as an initial design tool. 


EEE RSS SS SS ee ee ae epee 

FORM PAGE 1 
DD /2""..1473 | ) UNCLASSIFIED 
S/N 0101-807-6811 Security Classification 


75 ‘A-31408 


UNCLASSIFIED 


Security Classification 


v 
KEY WOROS 


PERTURBATIONS 

EARTH-MOON SYSTEM 

MATCHED ASYMPTOTIC EXPANSIONS 
TRAJECTORIES 


FoR 1473 (BACK) 


0t-807-6823 76 ncn Ass tate —______ poiFiiD A-31409 








Cc 





